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Abstract. I explore the observational appearance of the photosphere of an ultrarelativistic flow with internal 
dissipation of energy ( "dissipative" GRB model). As a case study, I use the magnetic reconnection model (AC 
model) that makes robust predictions on the energy dissipation rates at different radii in the flow. With analytical 
and numerical tools for the radiative transfer problem, I show that the flow develops a hot photosphere where 
inverse Compton scattering leads to highly non-thermal spectrum. For a wide range of luminosities and baryon 
loadings of the flow, this spectrum is very close to the observed prompt GRB emission. Its luminosity ranges from 
~ 3 to 20% of that of the total energy input. 
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1. Introduction 

Although much progress has been made the recent years 
in our understanding the central engine and of the after- 
glow emission of Gamma-ray bursts (GRBs), much less 
is known about the physical processes responsible for the 
prompt emission. Theoretical arguments related to the so 
called compactness problem (e.g. Piran 1999) suggest that 
the emitting material is ultrarelativistic with bulk Lorentz 
factors r > 100. The typical observed GRB spectrum has a 
characteristic peak (in E- f(E) representation) in the sub- 
MeV range that is smoothly connected to low and high fre- 
quency power-laws, usually modeled with the Band spec- 
trum (Band et al. 1993; Preece et al. 1998). Models for 
the prompt GRB phase have to account for both the ac- 
celeration of the flow and its observed spectral properties. 

For a flow to be accelerated to high bulk Lorentz fac- 
tors, it must start with high energy-to-rest-mass ratio. 
Depending on whether the energy is in thermal or mag- 
netic form, one has a fireball (Paczynski 1986; Goodmann 
1986) or a Poynting-flux dominated flow (Thompson 1994; 
Meszaros & Rees 1997; Spruit et al. 2001; Drenkhahn & 
Spruit 2002; Lyutikov & Blandford 2003). In the fireball 
model, the flow passes through an initial phase of rapid 
acceleration where T ~ r until most of the internal energy 
has been used to accelerate the flow (unless the baryon 
loading is very low and radiation decouples from matter 
before the acceleration phase is over). Further out, inter- 
nal shocks (Rees & Meszaros 1994; Sari & Piran 1997) can 
dissipate part of the kinetic energy and power the prompt 
emission. In magnetic models for GRBs, the flow accel- 
eration is more gradual (e.g. Drenkhahn & Spruit 2002; 
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Vlahakis & Konigl 2003) and dissipation of magnetic en- 
ergy through magnetic reconnection (Drenkhahn 2002) or 
current driven instabilities (Lyutikov & Blandford 2003; 
Giannios & Spruit 2006) can directly power the prompt 
emission. 

If the dissipated energy leads to fast particles, then 
the synchrotron or synchrotron self Compton mechanism 
appears as a natural one for the prompt emission. On 
the other hand, the large number of bursts with low en- 
ergy slopes steeper than the synchrotron model predic- 
tions (e.g. Crider et al. 1997; Frontera et al. 2000; Girlanda 
et al. 2003), have triggered alternative suggestions for 
the origin of the GRB emission. These include saturated 
Comptonization (Liang 1997), Comptonization by ther- 
mal electrons (Ghisellini & Celotti 1999) and the photo- 
spheric emission of the flow (Meszaros & Rees 2000; Ryde 
2004). 

Both fireball and magnetic models predict some degree 
of photosphcric emission that comes from the region where 
the flow becomes optically thin to Thomson scattering. If 
no energy dissipation takes place in the photospheric re- 
gion, this emission is expected to be quasi-thermal and 
its peak should be clearly observed in the case of fireball 
models (Daigne & Mochkovitch 2002); which is not the 
case. On the other hand, if there is substantial energy dis- 
sipation (through internal shocks, magnetic reconnection 
or MHD instabilities) close to the photospheric region, the 
emitted spectrum can strongly deviate from a black body 
(Meszaros & Rees 2000; Pe'er et al. 2005, 2006; Ramirez- 
Ruiz 2005). 

In this work, I study the spectrum of the photosphcric 
emission in dissipative GRB models. I use the magnetic re- 
connection model (Drenkhahn 2002; Drenkhahn & Spruit 
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2002) since it makes robust predictions for the characteris- 
tics of the flow and the rate of magnetic energy dissipation 
at different radii and study in detail the radiative trans- 
fer in this model. Similar considerations can be applied to 
other dissipative models (see also Pe'er et al. 2006). 

The study shows that while deep inside the flow ra- 
diation and matter are in approximate thermodynamic 
equilibrium, direct heating of the electrons through en- 
ergy dissipation leads to increase of their temperature as 
a function radius at radii where the flow is still thick to 
Compton scattering (about a factor of ten below the pho- 
tospheric radius). 

A Monte Carlo code has been implemented to study 
the effect of Compton scattering at the photospheric re- 
gion. The electron temperature is determined by bal- 
ancing the heating by dissipation with the cooling by 
Comptonization. The calculations show that the flow de- 
velops a hot photosphere with comoving electron temper- 
atures of the order of ~ 30 keV (and larger further out) for 
a large range of luminosities and baryon loadings of the 
flow. Inverse Compton scattering of the underlying ther- 
mal radiation leads to a Band-like spectrum which peaks 
at ~ 1 MeV (in the central engine frame) followed by a 
nearly flat high energy spectrum (with E ■ f(E) ~ E°) 
that extends up to a few hundred MeV. 

In the next Section, I discuss various dissipative mod- 
els for GRB outflows and review the main features of the 
magnetic rcconnection model. Analytical estimates for the 
radiative transfer and the electron temperature in the pho- 
tospheric region are given in Sect. 3. The description of the 
Monte Carlo code and the numerical results are presented 
in Sect. 4. The last two Sections contain the discussion 
and the conclusions. 

2. Dissipation close to the photosphere 

The fireball model for GRB outflows (Paczynski 1986; 
Goodman 1986) predicts the emission of a photospheric 
component from the region where radiation and mat- 
ter decouple. This decoupling takes place close to the 
Thomson photosphere of the flow. The photospheric ra- 
diation can carry away most of the luminosity of the flow 
(if the latter is still in its acceleration phase when the de- 
coupling takes place) or a smaller fraction <~ (r s /r p h) 2,/3 
of it depending on the baryon loading of the flow. Here 
r s stands for the saturation radius where the flow reaches 
its terminal Lorentz factor and r p h is the photospheric ra- 
dius. In the absence of energy dissipation around r p h, the 
photospheric component is expected to be quasi-thermal. 
Although a small fraction of GRBs do show thermal emis- 
sion (Ryde 2004) , the typical prompt GRB spectrum can- 
not be fitted by a black body. 

Internal shocks (Rees & Meszaros 1994) in an unsteady 
flow are usually invoked to dissipate a fraction of the bulk 
kinetic energy of the flow into fast particles and random 
magnetic fields and to power the prompt emission. For a 
range of the parameter space of the internal shock model, 
these shocks can take place below the Thomson photo- 



sphere, influencing the strength and the spectrum of the 
photospheric component (Rees & Meszaros 2005; Pe'er et 
al. 2005, 2006). 

In the case of Poynting-flux dominated flows, both ac- 
celeration of the flow and emission of radiation can take 
place through dissipation of magnetic energy. Magnetic 
energy can be dissipated directly through rcconnection in 
a flow where the magnetic field changes polarity on small 
scales (AC flow; Drenkhahn 2002; Drenkhahn & Spruit 
2002) or via MHD instabilities in an axisymmctric flow 
(Giannios & Spruit 2006). The magnetic dissipation and 
acceleration in these models is gradual and takes place 
over several decades in radius that typically include the 
photospheric radius. The photospheric component can be 
quite strong, of the order of <20% of the luminosity of the 
flow. 

Here, I focus on the observational appearance of the 
photospheric component in an "AC" flow. The basic fea- 
tures of the dynamics of the flow are reviewed in the 
next Section, while the detailed dynamical calculations are 
presented in Drenkhahn (2002) and Drenkhahn & Spruit 
(2002). 

2.1. The reconnection or "AC" model 

An important physical quantity of the flow in the "AC" 
model is the ratio cr of the Poynting flux to kinetic en- 
ergy flux at the Alfven point tq. This is a quantity that 
plays similar role to the "baryon loading" parameter that 
appears in fireball models, the ratio <r determines the 
terminal bulk Lorentz factor of the flow which is of 

3/2 

the order of <~ <j . The flow must start Poynting-flux 
dominated with cr > 30 for it to be accelerated to ultra- 
relativistic speeds with > 100. 

In an "AC" model, such as produced by an inclined 
rotator (Coroniti 1990), the magnetic field in the flow 
changes polarity on small scale A of the order of the light 
cylinder in the central engine frame (i.e. A ~ 2tvc/Q,, where 
n is the angular frequency of the rotator). In the flow, 
the energy density of the magnetic field is larger than 
the rest mass energy density and the reconnection speed 
is increased by the relativistic kinematics to subrelativis- 
tic speeds (Lyutikov & Uzdensky 2003; Lyubarsky 2005). 
Magnetic reconnection is modeled to proceed at a frac- 
tion e ~ 0.1 of the Alfven speed (which is essentially the 
speed of light for a magnetically dominated flow). Because 
of dissipation of the magnetic energy, about half of the 
Poynting flux converts into kinetic flux (i.e. acceleration 
of the flow) and the other half into internal energy of the 
flow. The dissipation is gradual and takes place up to the 
"saturation radius" r s where reconnection stops and the 
flow achieves its terminal Lorentz factor. 

Under the assumption of 1-D, steady flow, the rela- 
tivistic MHD equations can be solved analytically for radii 
ro <C r <C r s , yielding the following self similar scalings for 
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the comoving number density and the comoving magnetic 
field strength (Drcnkhahn 2002) 
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respectively. The notation A — \Q X A X is used and the 
"reference values" of the model parameters are Co = 100, 
e = 0.1, n = 10 4 rad s _1 , L = 10 52 erg s _1 sterad -1 ; very 
close to those used in previous studies of the model (see 
Drenkhahn 2002; Drenkhahn & Spruit 2002; Giannios & 
Spruit 2005). 

The model predicts gradual acceleration of the flow 
with r ~ r 1 / 3 in the regime ro C r < r s , while no further 
acceleration takes place above the saturation radius. The 
bulk Lorentz factor of the flow is, thus, approximately 
given by the expression 

1/3 
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The saturation radius is given by 
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Another characteristic radius of the flow is the 
Thomson photosphere. The optical depth between two 
points in the flow can be found by integrating the ex- 
pression dr = T(l — (3 cos 8)n' cttcIs (Abramowicz et al. 
1991), where /3 = v/c and 9 the angle of a photon path 
with respect to the radial direction. Integrating the previ- 
ous expression from r to oo for a radially moving photon, 
one gets the characteristic Thomson optical depth as a 
function of radius 
20 L 52 



2/3 5/2 • 
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The location of the photosphere is found by setting r = 1 
and solving for r 



r ph,n 
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The location of the photospheric radius with respect 
to the saturation one determines the nature of the emitted 
spectra expected from the flow. If, for example, r p h 3> r s , 
all the energy dissipation takes place in optically thick 
conditions and the radiation is efficiently thermalized. On 
the other hand, if r s > r p h, energy dissipation at moder- 
ate and low optical depths leads to a photospheric emis- 
sion that has a highly non-thermal appearance that can 
be directly related with the characteristic GRB emission. 
These points become clear in the next sections where the 
radiative transfer around the region of the photosphere is 
studied with analytical and numerical tools. 



3. Analytical estimates 

Energy dissipation in an "AC" flow is gradual and takes 
place over many decades of radii up to the saturation ra- 
dius r s . Half of the dissipated Poynting flux serves to ac- 
celerate the flow, while the other half is released as inter- 
nal energy in the flow. The rate of energy density release 
in a comoving frame can be found by the following con- 
siderations. The time scale over which the magnetic field 
decays is that of advection of magnetic field of opposite 
polarity to the reconnection area. The reconnection speed 
is v r = eva — ec, while the magnetic field changes po- 
larity over a length scale A' = 27iTc/(£n) measured in a 
frame comoving with the flow. The decay timescale for the 
magnetic field, therefore, is 

A' 2tiT 

idec — — — — pr- (') 
v r ell 

Using the last expression and eqs. © and the rate 
of magnetic energy density dissipation in the comoving 
frame is 

B' 2 /8tt 1.5 ■ 10 15 L 5 



^52 
3/2 
7 0,2 



ergem 3 s . 
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The fate of the released energy is rather uncertain. An 
interesting possibility is that dissipation leads to MHD 
turbulence where particle acceleration can take place by 
scattering of photons by Alfven waves at the r ~ 1 re- 
gion of the flow (Thompson 1994). On the other hand, 
the magnetic energy can directly be dissipated to the par- 
ticles in the flow, most likely to the electrons due to their 
higher mobility. Here I assume that a fraction f c of order 
unity of the energy is heating the electrons of the flow. 
The electrons are assumed to have a thermal distribution. 
Deep in the flow, the relaxation timescale of the electrons 
due to Coulomb collisions is fast enough to ensure ther- 
malization. I return to this issue in Sect. 4.3. 

Deep inside the flow the released energy is efficiently 
thermalized and shared between particles and radiation. 
Assuming complete thermalization, integration of the en- 
ergy released at different radii in the flow, taking into 
account adiabatic cooling, leads to the following expres- 
sion for the comoving temperature of the flow (Giannios 
& Spruit 2005) 

-1/4 

7k = , 1/9 keV. (9) 
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Comoving temperatures of a few hundred eV are typical 
at the region where the optical depth of the flow is of order 
of unity (see the last expression and eq. 10). At these tem- 
peratures the Compton scattering cross section dominates 
over free-free absorption so that deviation from thermal 
spectra is expected. It turns out, however, that the as- 
sumption of equilibrium of radiation and particles breaks 
down deeper in the flow and that the "photospheric" spec- 
tra are highly non-thermal for typical GRB parameters. 

The last point becomes clear if one looks at the en- 
ergy balance of the electrons in the flow. The electrons 
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are heated by magnetic energy dissipation and cool mainly 
through radiative cooling (the adiabatic cooling and the 
energy exchange with the protons are much slower). An 
obvious candidate for radiative cooling of the electrons in 
a strongly magnetized flow is synchrotron cooling since the 
energy density of the magnetic field is larger than that of 
radiation 



U B B' 2 /8ir 



Ur 



23 



(10) 



This suggests dominance of synchrotron over Compton 
cooling if the flow is optically thin to synchrotron emis- 
sion. It turns out that this is not the case in the photo- 
spheric region under consideration and that synchrotron 
radiation is strongly self-absorbed (I return to this issue 
in the end of this Section). The dominant cooling process 
for the electrons is, thus, Compton cooling. 

The Compton cooling rate for the electrons is given by 
the expression (e.g. Rybicki & Lightman 1979) 

ffcomp = 4n'e c CCr T C/r, (11) 

where O = i^sT e /(m e c 2 ) and ctt is the Thomson scat- 
tering cross section. Equating the heating rate (given by 
eq. ||HJ|) to the cooling rate of the electrons, I solve for the 
electron temperature 
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where / e> i stands for f e = 1. So the electron tempera- 
ture increases as a function of radius and equilibrium of 
the electrons and radiation is reached deep in the flow at 
the radius where T t h = T c . The last equation defines a 
characteristic radius which I call the equilibrium radius 



r eq,H 
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where eqs. (0 and i|12[l have been used to derive the last 
expression. Note that equilibrium of matter and radiation 
is achieved at a radius that is factor of ten shorter than 
that of the Thomson photosphere, or at an optical depth 
(combining eqs. © and (fHl ') 



46 
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which depends weakly on the model parameters and is 
much larger than unity. 

At radii r > r oq , the electron temperature is higher 
than that of the radiation field and upscattering of the 
photons takes place. Because of the increase of the elec- 
tron temperature with radius, this upscattering becomes 
more efficient close to the location of the photosphere. The 
electron temperature there can be found by combining eqs. 
JfjJ) and i|12|) which give Tc.ph = 40/ Cj i keV independently 
of the parameters of the model except to the fraction f e . 
For temperatures T c > 40 keV in the < 1 region, the 



characteristic signature of unsaturated Comptonization is 
expected to lead to spectra with hard non-thermal appear- 
ance. 

The previous estimates are based on the assumption 
that synchrotron emission in the flow is strongly self- 
absorbed and therefore the electrons cool mainly through 
Compton upscattering the radiation field. I can now check 
the validity of this assumption by estimating the syn- 
chrotron (and free-free) absorption optical depths in the 
flow and the characteristic turn-over frequency below 
which radiation becomes optically thick. The absorption 
optical depth for a photon traveling from radius r to 2r is 
given by the expression r v = a v r/T, where a„ is the ab- 
sorption coefficient as measured by the comoving observer. 
Note that I have not used the "effective" optical depth 
even though the inner parts of the flow are Thomson thick. 
This is justified by the relativistic nature of the problem. 
A photon that is emitted at a radius r that corresponds to 
an optical depth r> 1, does not undergo ~ r 2 scatterings 
until it escapes as it would be expected from random walk 
arguments. It suffers ~ r scatterings since it preferentially 
scatters along the radial direction (with an typical angle 
9 ~ 1/r <C 1 with respect to the radial direction) in the 
central engine frame. 

For a thermal plasma and for characteristic comoving 
photon energy hv <C feB^e, the absorption coefficient is 
related to the emission j v through the well known expres- 
sion a v = j v /2v 2 m c Q c . For the mildly relativistic plasma 
under consideration, the synchrotron emission is well ap- 
proximated by eq. (13) of Wardzinski & Zdziarski (2000; 
see also Petrosian 1981; Petrosian & McTiernan 1983). 
Using their expressions for synchrotron emission, I have 
calculated the turnover frequency defined by the ex- 
pression T Vt = 1 at different radii in the flow and have 
verified that v t ^> v c (where v c = eB 1 /2irm c c is the cy- 
clotron frequency) for a large range of the model parame- 
ters and up to radii that correspond to low optical depths 
~ 0.01 (or corresponding radii given by eq.[SJ); sufficient for 
this study. Below v t the synchrotron emission is strongly 
self-absorbed, resulting in synchrotron cooling rates in the 
plasma orders of magnitude less than the Compton cool- 
ing rates. 

As the analytic estimate Ijl2|l shows, the electron tem- 
perature increases rather fast and can upscatter the pho- 
ton field causing deviations from thermal spectra at a ra- 
dius of about a factor of ten below the Thomson photo- 
sphere. On the other hand, this estimate has its limitations 
since it gives the electron temperature by balancing the 
electron heating of the electrons to the Compton cooling 
rate assuming thermal distribution of photons. A more 
accurate calculation involves a self consistent determina- 
tion of the electron temperature simultaneously with the 
actual distribution of the photon field. This calculation is 
the topic of the next Section. 
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4. Numerical Study 



To study the emergent spectrum and the electron temper- 
ature above the "equilibrium" radius in detail, a Monte 
Carlo code has been developed. It simulates the Compton 
scattering in a flow with density and bulk Lorentz fac- 
tor given by eqs. and 10} respectively, by following 
the scattering random walk of a large number of photons 
(Pozdniakov et al. 1983). The special relativistic effects 
related to both the bulk motion and the scattering cross 
section of photons in the flow are taken into account. The 
inner boundary of the computational domain is taken at 
the last radius where the electrons are in equilibrium with 
the radiation field, i.e. the equilibrium radius defined by 
eq. (fT3")) . At this radius a black body photon spectrum 
is injected with comoving temperature given by eq. @ 
evaluated at r oq (given by eq. JT3J)) 



1000 F 



T 
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/eT 7 (^)3 /9 -o4 18 



-2/27 
J 52 



keV. 



(15) 



The outer boundary of the computational domain is at 
r ou t where the optical depth becomes small (the value r = 
0.1 is used). 

Obviously the result of the Comptonization depends 
critically on the electron temperature which, in turn, de- 
pends on the Compton cooling rate. The accurate deter- 
mination of T e (r) demands an iterative method where the 
flow is divided in a large number of spherical shells (of 
the order of 10 2 ) and an initial guess for T e (r) (e.g. the 
analytic estimate of eq. (|12|) ') is given. Then T c is var- 
ied in the different shells until balance of the heating and 
cooling rates is achieved at all radii. This method is, how- 
ever, rather computationally demanding. Fortunately, it 
turns out that one can get very accurate results assuming 
a power-law dependence of the electron temperature with 
radius of the form 



To -T 



(16) 



where T and s are two free parameters to iterate so that 
heating approximately balances cooling at every radius. 
The results of such iteration and the emitted photospheric 
spectra and their energetics for various values of the pa- 
rameters of the model are the subject of the next Section. 

4.1. Results 

The first step toward calculating the photospheric spectra 
is to determine the electron temperature parameters To 
and s for which Compton cooling balances the heating 
rate of the electrons at every radius in the flow. A first 
guess for To and s comes from the analytic estimate l|12|) . 

In figure 1, 1 plot the resulting electron temperature as 
a function of radius for different values of the parameters. 
At small radii, radiation and matter are in thermal equi- 
librium. Above the equilibrium radius introduced in eq. 
1|13(1 . the electron temperature increases following the ex- 
pression (|16|l for the numerically computed parameters Tq 




le+10 



le+ll 

radius r(cm) 



le+12 



le+13 



Fig. 1. The comoving electron temperature as a function 
of radius from the numerically determined energy balance 
with the radiation field. The baryon loading of the flow 
is given in the legend. The parameters that are not given 
are taken to have their reference values (see Sect. 2.1). 
For small radii radiation and the electrons are in thermal 
equilibrium. Triangles mark the radius where radiation 
and electron temperature get out of equilibrium, aster- 
isks mark the photospheric radius and circles mark the 
saturation radius. Note that the electron temperature at 
the photosphere is ~ 25 keV for all the parameter range, 
except from the high baryon loading (or low ctq) cases. 



and s which are close to the analytic estimate T oq and 5/3 
respectively. The asterisks in Fig. 1 show the location of 
the Thomson photosphere. Note that the temperature at 
the location of the photosphere is rather independent on 
the parameters of the model and about ~ 25 keV, slightly 
lower than the analytical estimate (|12fl . For the high <jq 
cases, the electron temperature keeps increasing until the 
outer boundary of the simulation. For moderate and high 
(To, the temperature is ~ 200 keV at this outer boundary. 
The very similar temperatures at the same optical depths 
in the flow lead to similar y-Comptonization parameter 
and, therefore, similar emitted spectra for a large range 
of the parameters of the model. This point becomes clear 
further down in this Section where the computed spectra 
are presented. 

Note that in the low do ( or large baryon loading) cases, 
there is a point above which the electron temperature 
drops again as a function of radius. This point is the satu- 
ration radius r s (defined by eq. iQJ) and marked by circles 
in Fig. 1). Above this radius no significant dissipation of 
magnetic energy takes place and the heating rate of the 
electrons drops quickly. The region r > r s does not con- 
tribute much to the emergent spectra (except from some 
degree of adiabatic cooling) since the electrons become 
essentially cold. 

In figure 2, I show the cooling and heating rate as a 
function of radius for the numerically iterated values of 
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Fig. 2. The heating and cooling rates, dashed and solid 
lines respectively, as a function of radius for different pa- 
rameters of the flow given in the legend. The rest of the 
parameters are assumed to have their reference values (see 
Sect. 2.1). The cooling and heating rates balance each 
other well, except from a narrow region close to the "equi- 
librium" radius in the inner boundary of the computation 
domain. This illustrates that a power-law is a good ap- 
proximation for the dependence of the electron tempera- 
ture with radius. 

To and s. Note that the power law modeling for the ra- 
dial dependence of the electron temperature (see eq. (|16[) ) 
leads to a very good balance except in the very inner re- 
gion of the computed domain. So the modeling of T c ap- 
pears to be accurate enough for the purpose of this calcula- 
tion. Obviously, a smoother fitting function could be used 
around the equilibrium radius, which would give a better 
description in this region. However, since the transition re- 
gion is rather narrow, only a small fraction of the energy is 
dissipated there. This simplification is not expected to in- 
troduce large errors to the calculations and more detailed 
modeling of the transition has not been pursued. 

The emerging spectra from our calculations are given 
in Fig. 3 in E ■ f(E) representation for different values of 
the parameters of the model. The spectra are plotted in 
the central engine frame with arbitrary normalization (the 
energetics are discussed in the next paragraphs). The peak 
of the spectrum is close to 1 McV for a large region of 
the parameter space which is relevant for GRB flows. The 
clustering of the peak around 1 MeV can be understood 
by the following considerations. The temperature of the 
photons at the equilibrium radius in the central engine 
frame is 

^^^(re^Te^lSO/^/ 9 ^)^ 3 ^/, 3 keV, (17) 

where in the last step the expressions (I15II . I@J and (|13[) 

have been used . The photon spectrum at r eq peaks at 

1 Note that this expression for the photon temperature does 
not apply for the very low o~o < 25 — 30 cases for which dissipa- 



Fig. 3. The resulting E ■ f(E) spectrum for different val- 
ues of the parameters of the flow (photon energies in the 
central engine frame). The moderate and high tro cases ex- 
hibit similar spectra with a peak that clusters at ~ 1 MeV 
and a flat high energy tail that is a result of unsaturated 
Comptonization near the photospheric radius. For low <jq 
dissipation stops below the photosphere and the resulting 
spectrum is quasi-thermal. 

- 4T™ - 600 keV (in a E ■ f(E) representation) for 
the reference values of the parameters and depends very 
weakly on the flow parameters (see eq. I|17|)l. The inverse 
Compton scattering that takes place above the equilib- 
rium radius leads to a moderate increase of the energy of 
the peak and, most important, to a power-law high energy 
emission with photon number index ~ — 2.3 that extends 
up to a few hundred MeV. This high energy part of the 
spectrum is the result of unsaturated Comptonization tak- 
ing place close to the Thomson photosphere and appears 
for flows with uq > 50. 

Both the peak and the high energy slope lie well within 
the observed range of the typical GRB spectrum, confirm- 
ing that the photospheric components may be responsible 
for the prompt GRB emission. The low energy slope of 
the spectrum depends on the energy range over which one 
attempts to fit it. Far from the peak, its slope is about 
f(E) oc E 2 (i.e. with photon number index a ~ 1). When 
fitted in the BATSE range (assuming a burst at z=2) with 
the Band spectrum (Band et al. 1993), typical values of 
a ~ —0.3 are found for the spectral slope below the peak 
frequency. These values are compatible with those mea- 
sured in hard bursts and cannot be explained with the 
simplest synchrotron models. 

In rather low gq < 50 flows, the emerging spectrum 
has much weaker emission above its peak. In these cases 
the magnetic dissipation stops close to, or even below, the 
Thomson photosphere and there is only weak Compton 
upscattering taking place in the photospheric region. 

tion does not proceed up to r cq which is the condition for the 
expressions @ and, therefore, 1151 and 11171 to be applicable. 
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Quasi-thermal emission has been observed in a fraction 
of GRBs (Ryde 2004; 2005) and they may be a result of 
low cr (i- e - high baryon loading). For even smaller values 
of <7 < 25 — 30, the magnetic dissipation stops bellow the 
equilibrium radius and the photon field suffers substantial 
adiabatic cooling before it decouples from matter. This re- 
sults in quasi-thermal photospheric spectrum that peaks 
at tens of keV or lower. 

The X-ray flashes are events that have spectral prop- 
erties very similar to these of the classical GRBs but with 
spectral peak below ~ 50 keV and are believed to be- 
long to the same family with GRBs (e.g. Barraud et al. 
2003). It is tempting to identify the X-ray flashes with 
these very low <jq flows. On the other hand, the spectrum 
of the X-ray flashes is similar to that of the classical GRBs 
which makes this identification unlikely. It appears more 
natural, in the context of this model, that X-ray flashes 
are typical GRB flows viewed off-axis (and therefore with 
low Speak)- Although the issue on the nature of the X-ray 
flashes is not settled, afterglow modeling seems to support 
this interpretation (e.g. Granot et al. 2005). 

4.2. Efficiency of the process 

A question quite relevant to the observational relevance 
of the photospheric component is its strength. A conve- 
nient quantity to measure the strength of this component 
is the photospheric efficiency e p h defined here as the ratio 
of the photospheric luminosity to the total luminosity of 
the flow. In Fig. 4, I plot the e p h for various <jq. For large 
(To, the efficiency is rather low <10% since most of the 
magnetic energy is dissipated further out in the flow. For 
moderate values of 00,2 ~ 0.5, the photospheric compo- 
nent becomes stronger with e p h ^ 15%, while its strength 
is reduced again for low <jq. In this case, the dissipation 
stops when the flow is still Thomson thick and radiation 
is cooled adiabatically before it decouples from matter. 
These numbers correspond to L 52 = 1, (eft) 3 = 1. Higher 
flow luminosities and values of efi lead to higher efficien- 
cies and vice versa (see Fig. 4). 

In Fig. 4, the "photospheric efficiency" in the BATSE 
range e p h,B is also plotted. The BATSE range is taken 
to be (30-2000) keV which correspond to the (30- 
2000) x (1+z) keV energy range in the central engine 
frame. For the redshift, I take z—2. In most of the param- 
eter space explored, more than ~60% of the photospheric 
emission is in the BATSE range; resulting in e p h,B close 
to e ph (see Fig. 4). 

In the high ao cases, additional radiation is expected to 
come from outside the computed domain since magnetic 
dissipation continues up to r s 3> r p ^. In this Thomson thin 
region, various of the assumptions underlying our calcula- 
tions may not hold anymore. I have assumed a thermal dis- 
tribution of the electrons which can no longer be justified 
outside the photosphere (see next Section). Furthermore, 
the electron temperature increases rapidly with radius. At 
some point it approaches that of the rest mass of elec- 




Fig. 4. The ratio of the photospheric luminosity in the 
whole energy range (black symbols) and in the BATSE 
energy range (red symbols) to the flow luminosity. Overall 
the photospheric radiative efficiency varies in the range a 
few to 20%. Higher flow luminosities and values of efl 
increase the radiative efficiency and vice versa. 



trons and pair creation should also be taken into account. 
Last, the non-thermal electron distributions are expected 
to lead to efficient synchrotron emission and synchrotron 
cooling cannot be neglected anymore. The problem of ra- 
diative transfer in this Thomson thin region depends crit- 
ically on the particle acceleration process in the regions 
where magnetic reconnection takes place. Some estimates 
on the processes at work in this region can be found in 
the Giannios & Spruit (2005). 

I discuss the issues of the thermalization of the electron 
distribution and of pair creation more quantitatively in the 
next subsections. 

4.3. Thermalization of the electrons 

In our study, I assumed a thermal distribution for the 
electrons. Particle acceleration in the current sheet in the 
reconnection regions can lead to highly non-thermal parti- 
cle distributions but, if the time scale of energy exchange 
through Coulomb collisions is shorter than the cooling 
timescale, thermalization of the electron distribution is 
achieved. The relaxation timescale for the electrons is 
given by (Spitzer 1967; Stepney 1983) 



_ So' _ in _ 8 29/6 JeA (. £il >3 g 0,2 



In A n'oTC 



,5/2 
J 52 



8,(18) 



where In A is the Coulomb logarithm. For the typical den- 
sities and temperatures of the problem at hand the value 
In A = 13 is used. I have also used the eq. and the 
estimate <|12|) for the electron temperature (the use of the 
more accurate numerically calculated values for T c leads 
to very similar results). 
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The cooling timescale of the electrons is 



t, 



^n'k B T c 



cool 



k 1n -7 7/3 (e^)3 7 V 

5 • 10 



0,2 



L 



r.2 



(19) 



-*Comp 

where the eqs. CD, © an d CJ are use d m tne 

last step. So, at radius of about 10 11 cm the relaxation 
timescale due to Coulomb collisions is shorter than the 
cooling timescale. However, the former has steeper radial 
dependence and the two timescales become equal at a ra- 
dius r* 



L 



r*ii 



3/5 
52 



(20) 



This radius is a factor of 2 inside the photospheric ra- 
dius for f e ~ 1 (compare the last expression with eq. ©)■ 
Thermalization of the electron distribution is therefore a 
fair approximation inside the photosphere. Above the lo- 
cation of the photosphere, deviations from thermal dis- 
tribution are expected. The actual electron distribution, 
however, is still expected (from energetic considerations) 
to peak at energies where heating and cooling rates bal- 
ance each other. 

4.4. Pair Creation 

In both the analytical estimates and the numerical inves- 
tigation, I have neglected pair creation. This is justified 
close to the equilibrium radius where the photon field is 
thermal with comoving temperature of about 1 keV; i.e. 
much below the pair creation threshold (see minimum of 
the curves in Fig. 1). However, the electron temperature 
increases further out in the flow, reaching values of the 
order of 200 keV near the end of our computational do- 
main at r = 0.1. A fraction of the photons are expected 
to be upscattered above the pair creation threshold in the 
comoving frame and lead to pair creation. To estimate 
how efficiently pair creation takes place, one can define 
the "comoving compactness" = n 7 <7Tr/r, where n 7 is 
the number density of photons that exceed the rest mass 
energy of the electrons (in the comoving frame) . It can be 
shown (Lightman 1982; Svensson 1982; Pe'er & Waxman 
2004) that if I' is much larger than unity, extensive pair 
creation takes place in the flow and the number of pairs 
can exceed that of the electrons related to baryons. On 
the other hand, if /' < 1 pair creation is not substantial. 

I have calculated n 7 directly from our Monte Carlo 
simulations at different radii and for different values of 
the parameters of the model. As expected, I' is found to 
increase with radius and to reach values up to I' < 0.2 for 
a large parameter space investigated: 0.3 < 00,2 — 
0.01 < L 52 < 10, 0.1 < (£fi) 3 < 10, < / e < 1. So the 
effect of pair creation is not expected to be important for 
parameters relevant for GRB outflows in the photospheric 
region under consideration. 

This result can also be understood in view of Fig. 3. 
The high energy part of the spectrum shows an exponen- 
tial cutoff at E c ~ 3 • 10 5 keV. In the comoving frame, 



this feature appears at an energy E c /T ~ m e c 2 . So only 
a small fraction (typically 10~ 3 ) of the photospheric lu- 
minosity is above the pair creation threshold, resulting in 
rather low compactness to pair creation. 

5. Discussion 

In this work, I have investigated the appearance of the 
photospheric component in Poynting-flux dominated GRB 
outflows in which the magnetic field changes polarity over 
small scales (AC model). An important characteristic of 
the model is that it predicts gradual dissipation of mag- 
netic energy through reconnection over many decades of 
radii including the region where the flow has Thomson op- 
tical depths of order unity (Drenkhahn 2002; Drenkhahn 
& Spruit 2002). 

Here I have shown that, if a large fraction of the en- 
ergy is dissipated directly to the electrons, the electron 
temperature increases rapidly at a wide region around the 
Thomson photosphere. The electron temperature is self- 
consistently calculated at different radii by balancing the 
heating and cooling rates. For this calculation both ana- 
lytical estimates and Monte Carlo simulations are used. 
It is shown that inverse Compton scattering of the un- 
derlying thermal radiation leads to spectra that peak in 
the 1 MeV range (in the central engine frame) and have 
power-law high energy part for a wide range of the model 
parameters. The high energy power law is a result is unsat- 
urated Comptonization that takes place at optical depths 
of order of unity. When fitted with the "Band" spectrum 
(Band et al. 1993), the spectra have low and high fre- 
quency spectral slopes and peak frequency in agreement 
with observations. This "photospheric" component is a 
significant fraction (from ~ 3% to more than 15%) of the 
luminosity of the flow and may, therefore, be responsible 
for the prompt GRB emission. 

For high baryon loadings (low do) in the flow, the en- 
ergy dissipation stops below the Thomson photosphere 
and no spectral component appears above the thermal 
peak. In this case, the photospheric emission is quasi- 
thermal and may be responsible for the appearance of a 
fraction of GRBs (Ryde 2004, 2005). 

Our calculation is limited to a region that extends up 
to a factor of ~ several above the photosphere (where 
the Thomson optical depth drops to r ~ 0.1). However, 
models with sufficiently low baryon loading (or equiva- 
lcntly high <7o) predict energy dissipation that continues 
further out in the flow. The spectra expected by these 
outer parts of the flow have been to some extent investi- 
gated by Giannios & Spruit (2005). On the other hand, 
as long as the issue of particle acceleration in the recon- 
nection regions is poorly understood, definite predictions 
on the relative importance of radiative mechanisms (e.g. 
Compton scattering, synchrotron emission) or the emitted 
spectra are hard to be made about the Thomson thin part 
of the flow. 

Our study was limited to the "AC" model because 
it makes robust predictions concerning the dynamics 
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of the flow and the rate of energy density dissipation. 
Furthermore, these predictions can take the form of simple 
algebraic expressions that simplify our study of radiative 
transfer in the flow. On the other hand, a large variety 
of models have or may have significant energy dissipation 
in the photospheric region that can lead to strong devi- 
ations from quasi-thermal spectra. In the context of the 
internal shock model, for example, shocks can also take 
place close to the location of the photosphere, leading to 
modifications of the photospheric emission (Meszaros & 
Rees 2000; Ryde 2004; Rees & Meszaros 2005; Pe'er et 
al. 2005). Another example of dissipative models is that 
of a strongly magnetized flow with an axisymmetric mag- 
netic field. Such flow is subject to current driven instabil- 
ities (kink instability) that lead to gradual dissipation of 
magnetic energy and rather strong photospheric emission 
(Giannios & Spruit 2006). 

6. Conclusions 

The standard fireball model for GRB flows predicts a 
rather strong photospheric component that is emitted 
when radiation and matter decouple. A photospheric 
component is also expected in magnetic GRB models 
(Drenkhahn & Spruit 2002; Giannios & Spruit 2006). In 
the absence of dissipative processes close to the photo- 
sphere, the photospheric component is expected to be 
quasi-thermal. In that thermal peak should be sys- 

tematically observed in the prompt GRB spectra, which is 
not the case (Daigne & Mochkovitch 2002). On the other 
hand, most of the GRB models for the prompt emission 
invoke energy dissipation at a large range of radii. If much 
energy dissipation takes place in the region of the photo- 
sphere, it can lead to large deviation of the photospheric 
component from purely thermal. Such investigation has 
been made in the context of internal shock and slow dis- 
sipation models by Pe'er et al. (2005). 

In this work, I have focused on the photospheric emis- 
sion expected from a strongly magnetized outflow in which 
the magnetic field changes polarity over small scales, re- 
connects and accelerates the flow gradually (Drenkhahn 
2002; Drenkhahn & Spruit 2002). Assuming that a large 
fraction of the dissipated magnetic energy heats the elec- 
trons in the flow, both analytical estimates and numerical 
calculations show that at a radius that lies a factor of ten 
below the Thomson photosphere r p h radiation and mat- 
ter are no-longer is thermal equilibrium. As a result of the 
energy dissipation, the flow develops a "hot" photosphere 
where electron temperature increases as a function of dis- 
tance with comoving temperatures ~ a few tens keV at 
the location of the photosphere. This result is rather inde- 
pendent of the model parameters such as the luminosity 
of the flow or the baryon loading ao- 

Inverse Compton scattering plays an important role at 
the photosphere, leading to Comptonization spectra that 
have characteristic non-thermal appearance. Fits of the 
numerically calculated spectra with the Band spectrum 
(Band et al. 1993) give parameters of the low/high fre- 



quency slope and the -E pea k of the spectrum in E-f(E) rep- 
resentation in agreement with observations. Furthermore, 
the observed clustering of the -E pea k in the sub-MeV range 
is a natural outcome of the model. The strength of the 
photospheric component is ~3-20% that of the luminos- 
ity of the flow and has most (>60%) of its energy in the 
BATSE energy range. I, therefore, conclude that the pho- 
tospheric component expected from magnetic models can 
to a large extent be responsible for the prompt GRB emis- 
sion. 

Acknowledgements. I thank Henk Spruit for many valuable 
suggestions and discussions and for carefully reading the 
manuscript. I acknowledge support from the EU FP5 Research 
Training Network "Gamma Ray Bursts: An Enigma and a 
Tool." 



References 

Abramowicz, M. A., Novikov, I. D., & Paczyhski, B. 1991, ApJ, 
369, 175 

Band, D., Matteson, J., Ford, L., et al. 1993, ApJ, 413, 281 
Barraud, C., Olive, J. -F., Lestrade, J. P. et al. 2003, A&A, 

400, 1021 
Coroniti, F. V. 1990, ApJ, 349, 538 

Crider, A., Liang, E. P., Smith, I. A., et al. 1997, ApJ, 479, 
L39 

Daigne, F., & Mochkovitch, R. 2002, MNRAS, 336, 1271 
Drenkhahn, G. 2002, A&A, 387, 714 
Drenkhahn, G., & Spruit, H. C. 2002, A&A, 391, 1141 
Frontera, F., Amati, L., Costa, E., et al. 2000, ApJS, 127, 59 
Ghirlanda, G., Celotti, A., & Ghisellini, G. 2003, A&A, 406, 
879 

Ghisellini, G., & Celotti, A. 1999, A&AS, 138, 527 
Giannios, D., & Spruit, H. C. 2005, A&A, 430, 1 
Giannios, D., & Spruit, H. C. 2006, A&A, 450, 887 
Goodman, J. 1986, ApJ, 308, L47 

Granot, J., Ramirez-Ruiz, E., & Perna, R. 2005, ApJ, 630, 1003 
Liang, E. P. 1997, ApJ, 491, L15 
Lightman, A. P. 1982, ApJ, 253, 842 
Lyubarsky, Y. E. 2005, MNRAS, 358, 113 

Lyutikov, M, & Uzdensky, D. 2003, Ap J, 589, 893 

Lyutikov, M., & Blandford R. D. 2003, astro-ph/0312347 

Meszaros, P., & Rees, M. J. 1997, ApJ, 482, L29 

Meszaros, P., & Rees, M. J. 2000, ApJ, 530, 292 

Paczyriski, B. 1986, ApJ, 308, L43 

Pe'er, A., & Waxman, E. 2004, ApJ, 613, 448 

Pe'er, A., Meszaros, P., & Rees, M. J. 2005, ApJ, 635, 476 

Pe'er, A., Meszaros, P., & Rees, M. J. 2006, ApJ, 642, 995 

Petrosian, V. 1981, ApJ, 251, 727 

Petrosian, V., & McTiernan, J. M. 1983, Phys. Fluids, 26, 3023 
Piran, T. 1999, Phys. Rep., 314, 575 

Pozdniakov, L. A., Sobol, I. M., Sunyaev, R. A. 1983, 

Astrophysics & Space Physics Rev. 2, 189 
Preece, R. D., Briggs, M. S., Mallozzi, R. S, et al. 1998, ApJ, 

506, L23 

Ramirez-Ruiz, E. 2005, MNRAS, 363, L61 

Rees, M. J., & Meszaros, P. 1994, ApJ, 430, L93 

Rees, M. J., & Meszaros, P. 2005, ApJ, 628, 847 

Rybicki, G. B. , & Lightman, A. P. 1979, .Radiative Processes 

in Astrophysics, Wiley, New York 
Ryde, F. 2004, ApJ, 614, 827 



10 



Dimitrios Giannios: Prompt emission spectra from the photosphere of a GRB 



Ryde, F. 2005, ApJ, 625, L95 
Sari, R., & Piran, T. 1997, MNRAS, 287, 110 
Spitzer, L. 1956, Physics of Fully Ionized Gases, Wiley, New 
York 

Spruit, H. C, Daigne, F., & Drenkhahn, G. 2001, A&A, 369, 
694 

Stepney, S. 1983, MNRAS, 202, 467 
Svensson, R. 1982, ApJ, 258, 335 
Thompson, C. 1994, MNRAS, 270, 480 
Vlahakis, N., & Konigl, A. 2003, ApJ, 596, 1104 
Wardzihski, G., & Zdziarski, A. A. 2000, MNRAS, 314, 183 



